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Single server queueing models are important in the study of a wide 
variety of stochastic systems. A particularly important example is the 
study of task schedules for computer systems with real-time applica- 
tions. In this paper, we present a class of approximations for the 
waiting time distribution in single server queueing systems with 
general independent (renewal) input and general (independent) ser- 
vice time distributions. These approximations allow the analyst to 
use as much (or as little) of the structure of the input and service 
processes as desired. Moreover, they also allow him to concentrate 
on specific quantities associated with the delay distribution that may 
represent the most relevant performance criteria. Examples include 
probability of delay, mean delay, and tails of the delay distribution. 
The methods given in this paper have been used to analyze the 
performance of task schedules in a variety of processor- based systems. 

I. INTRODUCTION 

Single server queueing models arise quite naturally in the study of 
a wide variety of stochastic server systems. A particularly important 
class is single processor computer systems, such as stored program 
control switching systems or nodes in a data communication network. 
In many such applications, it is important to keep as much of the 
structure of the interarrival and service time processes as possible in 
order to obtain realistic results; that is, the simplifying assumptions 
of exponential distributions cannot be made. While the resulting 
GI/G/1 queue may be extremely difficult to analyze, one is often 
content with reasonable approximations that incorporate the main 
features of the problem. In addition, one often desires results that are 
reasonably simple analytically, since the behavior of the GI/G/1 
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queueing model may be the input to the analysis of a more complex 
system. 

Our main purpose here is to present a class of approximations for 
the waiting time distribution, W(x), for such GI/G/1 systems which 
allows the analyst to use as much (or little) of the structure of 
the relevant input and service processes as desired. The resulting ap- 
proximations can be extremely simple in form (e.g., a single exponen- 
tial) or, with additional effort, more complex. In particular, we 
give relatively simple expressions for constants C and a, such that 
Wa(x) = 1 — Ce~ ax provides a good fit to the probability of delay, 
P D = [1 — W(0)] and the mean delay, w. Similar approximations are 
developed which provide a good fit to the tails of the delay distribution. 
Thus, these results can be used to study systems with a wide variety 
of performance criteria. Our results also lead to some new heavy traffic 
approximations, as well as a simple approximation for light traffic. 

While we feel that the approximations presented here are of use in 
themselves, we hope that the results obtained will stimulate more 
active research into the general method presented. Of particular inter- 
est is the problem of obtaining more quantitative error bounds to guide 
the user in the application of these techniques. 

We note that the methods given here have been used to analyze 
some rather complex GI/G/1 queueing systems which have arisen in 
the study of a certain class of computer systems with real-time appli- 
cations. 1,2 

Some key results are summarized in the next section to give a 
general idea of the nature of this work. Our basic approach is intro- 
duced in Section III where several (single) exponential approximations 
are derived. In Section IV, we look briefly at the implications of the 
approximations for heavy traffic. In Section V, we look at some special 
cases which lead to analytic statements about the accuracy of our 
various approximations, and in Section VI, their accuracy is assessed 
via several numerical examples. Extensions to more general functional 
forms are considered in Section VII and illustrated via another nu- 
merical example. Some final remarks are given in Section VIII. A 
summary of our notation, key formulas, and approximations are given 
in the Appendix for easy reference. 

II. SUMMARY OF SOME KEY RESULTS 

The main idea of this work is to assume a functional form for an 
approximation Wa(x) to the true waiting time distribution W{x) for a 
GI/G/1 queue. The Lindley integral equation for W(x) is then used to 
obtain the unknown coefficients. We obtain several approximations by 
assuming W(x) a W A (x) = 1 - Ce _<M: . A key approximation of this 
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type, referred to as approximation A * has C and a determined by eqs. 
(7) and (9) with a = a E (also see pages 299 and 300): 

c c * l - K{0) 

1 - K-{a E ) 

and 

C C* u + 



a a* K(0)-K-(a E ) 

where K(u) is the distribution function for the difference, w, + between 
the service time and interarrival time and 

K-(s) = e su dK(u) 



u+ = I udK(u) 
and a,E satisfies the characteristic equation 

e aEU dK(u) = 1. 



(The equation numbers of this section correspond to those of later 
sections.) Limiting properties and special cases for this approximation, 
as well as others are discussed. Several numerical examples are also 
given to illustrate the accuracy of these approximations. For an ex- 
ample H 2 /E VI system, the maximum relative error in the approxi- 
mation to Pd given by A* is found to be 5.0 percent, while the 
maximum relative error in the approximation to w given by A* is 
found to be 1.5 percent. Perhaps the simplest approximation consid- 
ered (using only the first two moments of u) is the one referred to as 
approximation A G ,o, which is also of the form Wa(x) = 1 - Ce~ ax , but 
where now the constants C and a are determined by 

-2w 

a= an = — —, (14) 

i.e., the heavy traffic a, and 

C = C G .o = [1 - Erf(|iZ|/>/2aK)]/[l + Erf(|u|/>/2aw)], (15) 

where Erf(*) is the error function. Note that (15) has the limiting 
behavior for small (u) (assuming o u tends to a finite limit) 



* The symbol * will be used to denote a random variable (see the Appendix for a 
summary of notations used). 
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C G ,o ~ l-2-v/-i^i, (16) 

U-.0 y 7T a u 

(p->l) 

which then results in the following limit for the approximate mean 
delay, Wgq 

Cg,o i/2 M - 

CLh u—0 CLh V 7T 

(p— 1) 

the right-hand side of eq. (17) is, in fact, a well-known lower bound for 
w (as u — > 0, e.g., see Ref. 4). 

The accuracy of these approximations, as well as several others we 
develop, are studied in considerably more detail in what follows. We 
also show how the approach taken leads to a class of approximations 
which can be made increasingly more accurate (with additional effort). 

III. DEVELOPMENT OF THE EXPONENTIAL APPROXIMATION 

We wish to obtain an approximate solution for the waiting time 
distribution, W(x), for a GI/G/1 queueing system. More specifically, 
we have a single server queue where the interarrival times, t, between 
customers are independent and identically distributed, taken from a 
general distribution, A(t) and the service times, t, are independently 
drawn from an arbitrary distribution B(t). A first in-first out discipline 
is assumed. If we let u = f — i and denote the distribution of u by 
K(u), then if an (equilibrium) waiting time distribution exists/ it 
satisfies the well-known Lindley integral equation 3 

W(x) =( W(x - y)dK(y), x > (1) 

(W(x) is equal to for x < 0). 

In this section, we consider approximations Wa(x) to W(x) of the 
form 

W A (x) = 1 - Ce~ ax . (2) 

Our objective is to determine suitable values for the constants C and 
a. For this purpose, we will use the relation (1). 

3. 1 Pointwise matching 

As is well known (see Ref. 5, pages 376 and 410), if the equation 

e au dK(u) = 1, (3) 



t Throughout we assume stability, i.e., that (1/a) = t < t = (1/X); p = (\/a) < 1. 
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possesses a real nonzero root, as, then 

1 - W(x) ~ Cc - "**. (4) 

Throughout this paper, we will tacitly assume that for all systems 
considered, eq. (3) does, in fact, possess such a root. Thus, if an 
exponential form is to be used, a£ is a reasonable choice for the 
parameter a in (2). On the other hand, substitution of (2) into (1) 
yields the relation 



1 - Ce~ ax = K(x) - Ce~ ax e ay dK(y) . (5) 

Hence, we see that if we choose a = as, (5) will be valid to order 
0(e~" EX )* as x tends to infinity. We can now determine C by requiring 
that (5) be valid at another value of x. For x such that [1 — K(x)] > 
we can solve (5) for C to obtain 

CM- [1 -* ( * )] *" ■ (6) 



1 - e ay dK(y) 

Now if we are interested in approximating the behavior of W(x) near 
x = 0, a reasonable choice might be to solve for C from (6) with x = 0; 
that is, choose 

c = Co(a)= i-«(°> = i-f(o> , (7) 

1 - K-(a) 1 - k-(a) 
where 

K-(s)= e su dK(u). 



[Note that clearly < C n (a) < 1.] 

We will refer to the exponential approximation of W(x) obtained by 
choosing a to satisfy (3) ["(5) at x = oo"] and C to satisfy (7) [(5) at 
x = 0] as approximation A E ,o, i.e., W Aeo = 1 - Co(a,E)e~ aEX . (The sub- 
script E refers to the use of the exact dominant root — we will consider 
other alternatives shortly.) 

Note that if the true delay distribution were given by a simple 
exponential, i.e., W(x) = 1 - Ce~ ax , then the value of C(x) determined 
from (6) would be identically a constant and W Aeq (x) = W(x). Thus, 
the variation of C(x) with x, as determined by (6), gives some indica- 



f That is, the difference between the left-hand side and right-hand side will tend to 
zero even when multiplied by Ce"'\ 
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tion of how exponential the waiting time delay distribution is. A more 
"sophisticated" approximation might be obtained by minimizing the 
variation of C(x) with respect to some suitable norm, however, we will 
concentrate here mostly on simpler techniques. 

Note also that C(x) as given by (6) appears in various (exponential) 
bounds for W(x). For example, Kingman 6 and Rossberg and Siegel 7 
prove the following inequality for the complementary waiting time 
distribution, W c (x) = 1 - W(x). 

C L e~ aEX < W c (x) < e~ aEX , 

where a E is given by (3) and 

C L = inf C(x). 

[1-K{x)]>0 



Ross 8 gives the bound 
where 



W c (x) < Cue- aEX , 



Cu = sup C(x). 

[1-KUI]>0 



(See also Ref. 9.) Hence, it seems quite natural to choose a specific 
value for x in C{x) to obtain an exponential approximation to W(x), 
with C L <C< Cu. 

3.2 An alternate method for determining C 

As we shall see shortly, A E ,o does provide a reasonable approximation 
to P D , but the resulting approximation for w is not as good. One might 
consider choosing C from (6) with another choice of x; however, here 
we consider another alternative. If we compute the mean waiting time, 
id, using eq. (1) we obtain 



w = xdW(x) = x 



W(0)dK(x) + d x W(x - y)dK(y) 



(8) 

Again, assuming the exponential form (2) for W A (x) [and, hence, 
that w A = (C/a)], (8) yields 

- = C - a+ - K'+(Q) (9) 

WA " a K(0) - K-(a) K-(0) - K-(a) * 

where u + = Jo udK(u), K + (s) = Jo e su dK(u) and £-(*) is as in (7). If a 
is chosen from (3) and C then from (9), we will refer to the resulting 
approximation asA E ,i, i.e., Wa e1 {x) = 1 - d(a E )e~ aEX . As we shall see, 
this will result in a good approximation for the mean delay via (9); 
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however, it is not clear from this relation that the resulting C need be 
less than 1. (It is clearly >0.) A correct probabilistic exponential form 
for Wa(x) which incorporates both (7) (and, hence, results in a reason- 
able approximation for the probability of delay) and (9) (and hence 
results in a good approximation for the mean delay) can be obtained 
in the following manner. First, find as from (3), then C = C* = Co(oe) 
from (7), and then the approximate mean delay, w — Ci{a.E)/a.E from 
(9). The constant a is then chosen to be 

. Co(a E ) 
a = a = r ( n \ aE - 

We will refer to this approximation as A*. We now look at a more 
organized way of obtaining such approximations. 

3.3 Moment matching 

Equation (9) was obtained by matching the first moment of the 
Lindley integral equation, i.e., from (8). We note that if we match the 
Oth moment in a similar way we are, in fact, also led to (7). That is, if 
we seek a solution to (1) of the form 1 — Ce -£ur , and match the first two 
moments (Oth and 1st) of (1) to determine both unknown constants, 
we are led to eqs. (7) and (9), i.e., 

..i-«w 



1 - K-(a) 

C _ u+ 

a ~ K(0) - it-(a) 

as a pair of equations in the two unknowns C and a. These can be 
combined and rewritten as 

= [K(0) - K-(a)][l - K(0)] 

u + [l - K-(a)] 

„ 1 - K(0) 

C = 7 r-^-. (11) 

1 - K-(a) 

Thus, one can first determine a from (10) and then C from (11). 
Clearly, if such an a exists, a E [0, oo], C E [0, 1]. 

Denoting the right-hand side of (10) by f(a), we note that 

f(a) -» 

o—n 

K(0)[l - KM] 
f(a) -> < oo 

and also 

QUEUEING SYSTEM APPROXIMATIONS 301 



-K'-(a)[l - K(0)f 

f (o) = ^ o — 

oT[l - K-(a)f 

has the properties 

f'(a) > 0, a > 

/'(a) -> 0. 

a—* en 

Moreover, 

rr.j " [1 " K(0)f ( ^- (a)[1 " *' (a)] + 2( ^ (a))2 < a>0 

' " + [1 - £-(a)] 3 

Thus, we see that (10) possesses a unique positive root, which can be 
obtained via the iteration scheme 

a n+ i = f(a n ) 

starting from any positive ao. 

Thus, we can always determine the desired root cto.i from (10) and 
then the resulting C = C 0>1 from (11). This approximation 

W A Jx) = 1 - C ,ie- aoii 

is referred to as Ao.i . 

3.4 Light traffic 

Under quite general conditions, the function K-(a) defined in (7) 
and with a satisfying (3) has the limit 

K-(a) -+ 0, 

e.g., fix B(t) and the "shape" of A(t) and let 1/A -► oo. Hence, for p 
small, (7) and (9) lead to the following intuitively appealing "light" 
traffic approximations {Alt) 

C lt = 1 - K(0) (12) 

^ = S = ^oT- (13) 

IV. HEAVY TRAFFIC 

In the case of heavy traffic, an obvious possibility for the exponent 
a is the heavy traffic limit (e.g., see Ref. 4): 

a„ = ^, (14) 

On 
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where u, at are the mean and variance, respectively, of u = t — t 
(service minus interarrival time). We can use this value for a in place 
of as in the approximations of Section III. In particular, if we use 
a = an in (7) we will refer to the resulting approximation as Ah.o- Note 
that the resulting C is strictly less than one unless o« = (p = 1). If we 
use a = an in (9), we refer to the resulting approximation as Ah,\. 
These two approximations (as well as those we discuss next) offer 
potential improvements to the more standard heavy traffic approxi- 
mation that can extend its applicability to lower load levels. 

In many cases, where heavy traffic approximations are applied, the 
structure of K{u) is either not known or too complex to be used in 
further analysis. Standard heavy traffic approximations generally only 
make use of the first two moments of K(u), i.e., u and at- We can also 
obtain results which use only u, at by making the standard heavy 
traffic assumption that K{u) is approximately Gaussian. For example, 
if we not only assume that a = an in (7), but further that K(u) is 
Gaussian, we obtain the following approximation for C using (7): 



l-Erff-"' 7 



Cg.o — • 



W£ 



O u 



1 + Erf 



{—) 

\vW 



(15) 



We refer to this approximation as A G ,o- Note that (15) implies the 
limiting behavior for | u \ small 



1- /-J*I 

'"-° 2 \u\ "—° ' w °" 

<p-i> 1 + \ \-- — -<p— i) 
V tt a u 



(16) 



Now combining (16) with (14), we obtain the following (limiting) 
approximation to the mean delay for | u \ small 



Cg.o 



1 2 

w G ,o = - ~ a "\/ _ - 

an u—o a« V it 

(p-«i) 



(17) 



In a similar manner, (9) yields the approximation Ag,\ 



w g .\ = 






2 - 



1 - Erf( ^- 



a a 



Erff-"' 7 



\y/2o u 



(18) 



which has the limit 
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u, Gi =^ _ JZL-ZlJ* ~ ±-±\E. (19) 

0,1 a// s^o 2|u| 2 \2 suo a w 2 V 2 

<p-»i) (p— i) 

Now if K(u) were, indeed, Gaussian, then the resulting (heavy 
traffic) mean delay, w G , is known to have the limiting bounds (e.g., see 
Ref.4). 

— -a u -v/-<«; G < — . (20) 

an V tr a.H 

Compare (17), (19), and (20). 

V. SOME SPECIAL CASES 

In this section, we look at the general behavior of these approxi- 
mations for some special cases where exact results are readily obtain- 
able. We begin by looking at the GI/M/l case where the approxima- 
tions of Section III are exact. This will serve as a simple illustrative 
example, as well as provide some additional insight into the behavior 
of the heavy traffic approximations of Section IV. 

5.1 Case:GI/M/1 
For a GI/M/l system, we can write (for u > 0) 

K(u) = | [1 - e- alu+t) ]dA(t) = 1 - pe~ au , u > 0, (21) 
Jo 

where l/a is the mean of the (exponential) service time and 



p= e-"'dA(t) - A(a). 

to 



u + = ccp xe- ax dx=p/a (22) 



We, thus, have 



K(0) = l-p. (23) 

Now if a satisfies (3), 



(«)- J 



K-(a) = e au dK(u) = 1 - e au dK(u). 



Hence, 



K-(a) = l-ap I e au e- au du = 1 4- -^— . (24) 

a — a 
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Using (23), (24) in (7) yields 



Co(a)=- — -, (25) 

a 



while (22), (23), (24) in (9) yield 

Ci(a) a — a 



aa 



(26) 



If a is chosen to be the exact root of (3), then these both lead to the 
exact waiting time distribution. (Recall, for a GI/M/1 system 

a 

e.g., see Ref. 10.) 

If we use the heavy traffic approximation for a, i.e., (14) then (25) or 
(26) can be used to obtain the following approximation for C 

0,0=^=1-^^ (27) 

a a o u p 

that is, a possible improvement to C = 1 for heavy traffic. Thus, the 
resulting mean delay is found to be 

, T , _ C *fl _ 1 1 _ a l l toa\ 

Wh.o Trprr . (28) 

ciH clh a 2\u\ a 

Note that this correction term to 1/a// as an approximation to the 
mean delay is typical. For example, for an M/G/l system it is easy to 
show that the true mean delay, w, satisfies 

- 1 l (! A 

W = T + - • 

a H 2 VA a 



5.2 Case:M/G/1 

For this case, we have 



K(u) = e M "- T) dB(T) = qe Xu , u<0, (29) 

Jo 

where 1/A is the mean of the (exponential) interarrival time and 



We, thus, have 



q = I e" AT dB(r) = B(\). 
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and 



u + = udK(u) = udK(u) - udK(u) 

JO J-oa •'-» 

-■-*£.-*■*- (s-M) (30) 

#(0) - q (3D 

Jt(a) = ?A e au e Xu du = -^-r • (32) 

J_ a + A 

Using (7), we obtain the following approximation (A e ,q) for the 
probability of delay, Pd 

C (a E )= — 1 ~~\ , (33) 

I—A" 

a^ + A 

while from (9) we obtain the following approximation (Ae.i) to the 
mean delay, w 

Ci(a E ) (A - a + qa)(a E + A) , Q . X 

We,i = = r • W 

as qaAa.E 

(Since we know that for an M/G/l queue, the true value of Pd is p, we 
could use this fact, together with (33) to obtain a simple approximation 
for a E ; we discuss this possibility shortly.) 

Using the known limits for the true a = a E as p — » 1, (33) readily 
yields the following: 

C (a E ) ~ p(l-p?a 2 ) ~ p (35) 

X-,0 X— 

(P-.0) (p—0) 

C (a*) ~ 1-—^—- -+ 1, (36) 

A— a 1 — q a A—a 

<p-.l> (p—1) 

where we have used the fact that for A -* 

f 00 A A 2 t^ 

o= e" AT dB(T) - 1-- + -T-. (37) 

Jo A -° a 2 

We similarly obtain from (34) 

_ CMe) A? / X\ A? 

"* — ^ x: 20^7) I w ^o 20^) (88) 

(p—0) (p— 0) 
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and for A — » a, clearly 



w E ,\ = 



Ci(a E ) 



ClE A—a ClE 

(p— 1) 



(39) 



Equation (38) shows that as A — * 0, we obtain the familiar P — K 
(Pollecek-Khinchin) formula for w in an M/G/l system, while (39) 
shows that with 

1 1 

a = a,E, we,\ ~ — ~ — . 
p—i as p-i oh 

Now recall for the approximation A*, we take C* = Co (as) and 

Co(a E )a E 



a T = 



Ci(a E ) 



Hence, for A*, the resulting exponent for the exponential approxima- 
tion is given by 

qaXaEil — q) 



a' = 



(A - a + qa)(a E + \ - q\) 



From (40) we obtain 



(40) 



(41) 



(42) 



where R is mean forward recurrence time of the service time distri- 
bution. Note we have used the fact that 

a E - 0. 

p— i 

The first of these is a very interesting relation. It is easy to show 
that for an M/G/l system, that as p — > the true dominant root a.E 
satisfies 



a* - 

A 
(p 


2 _ 1 

-0) 




a* - a E \ 1 - 

A-« [ 


a E (l - q) 
a 


\ - a E , 

J A-~„ 


<p— i) 




(p— 1) 



B(-a E ) = 



7 dB(r) 



(43) 



Hence, for example, if B(s) is rational, as tends to the smallest (in 
magnitude) pole of B(s), which is not consistent with (41). However, 
it is easy to show that for an M/G/l system, the mean delay condi- 
tioned on being delayed, is just the mean forward recurrence time, R — 
the dominant contribution is from the case where the arrival finds one 
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customer in the system. That is, for p -> 0, the dominant root really 
may not be the best quantity to use for an approximation, indeed, the 
limit (41) is most likely preferable. We, thus, have via A* an exponen- 
tial approximation to the true delay distribution for p -> which can 
be decidedly nonexponential but with the correct limiting probability 
of delay and mean delay. 

The approximation Aoj also has some very interesting properties for 
the M/G/l case. Considering (7) and (9) as denning two equations for 
the two unknowns C and a, we obtain the two relations 



and 



(a + A)= ^_. (45) 

A — a + qa 



Combining (44) and (45) we find that 



C = - = P . (46) 

a 



That is, for the M/G/l case, A ,i yields an exact expression for the 
probability of delay, P D . Using (44) or (45), we obtain the following 
approximation for the dominant root: 

ap 2 q pad - q)d - p) u ~ 

— ap = — — . (4/) 



a (p-(l-q)) "* p-d-q) 

Note that 

a(l-q)(l-p) 
a 

P-i ' Q 

but 



(48) 



2 



♦0 „ T 2 



as with A*. 

Note that since for an M/G/l system, we know that P D = p, we 
could have used this fact directly in either of the two approximations 
[(33), (34)] to obtain an approximation for the true dominant root, a E . 
This, of course, would result in 

pad-qHl-p) 

aE ri ^ — ' 

p-d-q) 

i.e., just (47). 
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For the mean delay, w, A .i yields the approximation 

p-(l-q) 







'■' «<1- 


q)d- 


-P)' 


which has the 


asymptotic properties 






Wo,\ 




Ar 2 






A— / 
(P-O) 21 1 ■ 


-P- 


n 




u>o.i 




Q 






X^a a(l - 

(p—1) 


-9)d 


-P)' 



(49) 



(50) 



(51) 



Comparing the asymptotic behavior of A 0A and A*, we see that while 
the former has a more desirable property for the approximation to P D , 
the behavior of the approximations to w and as are more desirable for 
A* than for A ,i. This leads us to consider another approximation, one 
that combines the best properties of each. The simplest method is to 
define approximation A** via wa-{x) = 1 — C**e~ a ** x , where 

C = Co.i 
and 

Co.i Co.iCtE 

a** = — — = — . 

w E ,\ Ce,\ 

The use of the heavy traffic oh, of course, yields similar results to 
those for Ae.o, Aba, ana< A* as p — » 1. Again, there is no simplification 
for the Gaussian case. 

VI. SOME NUMERICAL EXAMPLES 

We consider several numerical examples of varying complexity to 
illustrate the accuracy of the various approximations. Because of the 
large number of possible combinations of approximations and quanti- 
ties of interest, we will not cover all possibilities for every example. 

6.1 Case:M/D/1 

For the M/D/l case, 

KM = m = ? - ts u) = j^L, »>?. 

Hence, we can readily obtain eq. (3) for the dominant root, a E 



/' 



e""dK(u)=— ^— e" T "=l. (52) 

a + \ 
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Rewriting this equation as 

\T(e a ~- 1) = ar, (53) 

it is clear that in addition to a = 0, there is always one positive real 
root, a E , for any finite X. Moreover, 

a E -> 0, cle -* °°. 

Xf= p _l A-.0 

(p—0) 

The other quantities needed for the various approximations are 

q = B(X) = e~ p (54) 

K{0) = q = e- p (55) 

qX e~ p X , CC v 

K-(aE)=- !i - T = —-T < 56 > 

oe + A as + A 

•-e-M)-('-^¥) » 
-e-9-'-{ 

a^ = i (59) 



and, hence, 



\ 



aw = -^ = 2 P (l_-p) j (60) 



With a £ from (53) and (54) to (60), we can compute the various 
quantities predicted by our approximations. For simplicity, we will 
assume t = 1 in the following so that X = p. 

6.1.1 Probability of delay, P D 

Table I shows the value of P D as predicted by the various approxi- 
mations. Recall that P D from Ao.i is exact for this case. We see, 
however, that Ae,q{A* ) produces a good approximation to this quantity. 
The value of P D predicted by A E ,\ is seen to be significantly poorer, 
particularly for p small. Actually, our heavy traffic approximation A H ,o 
does significantly better than Ae,\ . Note that A G .o provides a reasonable 

Table I— Probability of Delay, P D , case M/D/1 



p 


True 


AoAA") 


A E .o(A*) 


Ah,0 


Ac.o 


Ae.i 


Alt 


0.01 


0.0100 


0.0100 


0.0100 


0.0149 


0.1920 


0.0327 


0.0100 


0.1 


0.1000 


0.1000 


0.0975 


0.1406 


0.2256 


0.1986 


0.0952 


0.5 


0.5000 


0.5000 


0.4756 


0.5647 


0.4462 


0.6170 


0.3935 


0.9 


0.9000 


0.9000 


0.8864 


0.8975 


0.8524 


0.9276 


0.5934 


0.99 


0.9900 


0.9900 


0.9884 


0.9885 


0.9842 


0.9928 


0.6284 
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Table II — Mean delay, w, case M/D/1 







A KA (A\ 














p 


True 


A**) 


■"O.l 


Ae.o 


A a j 


Ac., 


Aht 


Alt 


0.01 


0.00505 


0.00505 


0.00507 


0.00154 


0.00759 


12.53 


50.51 


0.00503 


0.1 


0.0556 


0.0550 


0.0565 


0.0270 


0.0832 


1.589 


5.556 


0.0536 


0.5 


0.5000 


0.4911 


0.5415 


0.3785 


0.7026 


1.033 


2.000 


0.3513 


0.9 


4.500 


4.478 


5.166 


4.279 


5.027 


4.895 


5.556 


0.8379 


0.99 


49.50 


49.56 


57.55 


49.35 


50.13 


49.88 


50.51 


0.9829 



approximation in the heavy traffic region and that A LT provides a good 
approximation in the light traffic region. 

6.1.2 Mean delays, w 

Table II shows the mean delays resulting from the various approxi- 
mations. We see that A E ,AA*) provides the best approximation to the 
true iv. The approximations A H ,\ and A .i are somewhat similar, with 
Ah,\ being better at higher loads, while A ,i appears better in the 
midrange. The approximation A E ,o appears quite poor for low loads 
but improves with increasing load. While not comparable with A H .i, 
Ag,\ does provide some improvement over A H t. In the light load region, 
Alt is seen to provide a good approximation. It is interesting to note 
that while A H ,u Ao,i, and Aht all use the heavy traffic a H from (60), 
which has the poor behavior 

a//-* 0, 

P-.0 

this only results in the anomalous behavior 

w A -> co 

for Ag.i and Aht- 

6.1.3 Dominant root, a L 

Table III shows the exponent that would be used in the (exponential) 
approximation for W(x) corresponding to each of the approximations. 
Recall that for Ae.o and Ae.i the exponent is exact, i.e., a = oe. We see 
that A* is somewhat better than A ,i for heavier loads and both tend 
to 2 as p — » — the inverse of the mean forward recurrence time of the 
service time. As should be expected, the heavy traffic a, a.H, is quite 
poor for light traffic. 

Table III — Dominant root, a, M/D/1 



p 


True 


A* 


A .t 


An 


0.01 


6.475 


1.974 


1.974 


0.0198 


0.1 


3.615 


1.775 


1.770 


0.1800 


0.5 


1.256 


0.9685 


0.9234 


0.5000 


0.9 


0.2071 


0.1979 


0.1742 


0.1800 


0.99 


0.0200 


0.0199 


0.0172 


0.0198 
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Note that the differences between P D and w as predicted by A E ,o 
and As.it particularly for light loads, indicates that the true delay 
distribution is essentially nonexponential. Yet, we see that assuming 
an exponential form can result in excellent approximations for P D and 
w. Moreover, the reasonable results predicted by Aha and Ah,\ for light 
loads where a H is far from any reasonable choice of a dominant root 
indicates that these techniques are somewhat robust with respect to 
the choice of the exponent, a. 

6.2 Case:M/E 2 /1 

We now look at several of the resulting approximations for the 
waiting time distribution, W A (x), and compare these with the exact 
results for an M/E 2 /l system. Specifically, we assume A(t) = 1 - e~ M 
and that B(t) is the convolution of two exponentials with unit mean. 
Hence, we have for the respective Laplace Stieltes transforms 

A + s 

1 

(1 + s)' 

For this system, the exact delay distribution, W(x) can readily be 
obtained via standard technique (e.g., see Ref. 10). The result is 

W(x) = 1 - Cie- a " - C 2 e- a * x , (62) 



B(s) = „ , . 2 . (61) 



where 



(2 - A) - [A(4 + A)] 
ai = 



1/2 



a 2 = 



2 
(2 -A) + [A(4 + A)] ,/2 



2 

(l-p)d-oi) 2 
Ci — 



C 2 = 



ai(a2 — cti) 
-(l-p)(l-a 2 ) 2 



ai(a 2 — cti) 
From (61) and (62), we can readily compute all of the quantities we 



need. We have 






q = &(\) = 


1 . ,_!_o. 


1 
u — 2 


(1 + A) 2 ' a "' 


A' 


and 


o2-2 + p, 





yielding expressions analogous to (54) to (60). 
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Tables IV, V, and VI show comparisons of the exact delay distribu- 
tion, as well as the mean delay, with some of the approximations we 
developed for p = At = 2A = 0.1, 0.5, and 0.9. The approximations for 
the probability of delay, Pd, and mean delay, w, show similar behavior 
to that for the M/D/l case. In addition, we see that the approximation 
Ae,i seems to have the best tail behavior (followed closely by A**). 

Thus, from these two examples, we might conclude that A** is the 
best overall approximation, while Ab,\ has slightly better tail behavior. 

6.3 Case: D/D 2 /1 (D/G/1) 

We briefly consider here the D/G/1 case and a simple numerical 
example which illustrates some interesting properties of our approxi- 
mations. 

For the D/G/1 case, we have 

,., ,B(u + d); u>-d 

""' 0; u<-d> 

Table IV — P(Delay < x) for case M/E 2 /1 



X 


True 


Ae.o 


A E .\ 


A ,i 


A' 


A** 




0.5 

1.0 

5.0 

10.0 


0.9000 
0.9220 
0.9414 
0.9962 
0.9999 


0.9014 
0.9323 
0.9534 
0.9977 
0.9999 


0.8760 
0.9148 
0.9414 
0.9971 
0.9999 


0.9000 
0.9257 
0.9448 
0.9950 
0.9997 


0.9014 
0.9268 
0.9460 
0.9950 

0.9997 


0.9000 
0.9261 
0.9454 
0.9951 
0.9998 


w 


0.1667 


0.1314 


0.1654 


0.1680 


0.1654 


0.1654 



Note: P = 0.1, True W(x) = 1 - 0.1667e-° 75001 + 0.06667e-' 200x . 

Table V — P(Delay < x) for case M/E 2 /1 



.r 


True 


A E .o 


Aba 


Am 


A* 


A" 




0.5 

1.0 

5.0 

10.0 


0.5000 
0.5644 
0.6272 
0.9084 

0.9848 


0.5119 
0.5922 
0.6593 
0.9192 
0.9866 


0.4666 
0.5547 
0.6277 
0.9117 
0.9854 


0.5000 
0.5742 
0.6374 
0.8998 
0.9799 


0.5119 
0.5859 
0.6488 
0.9058 
0.9818 


0.5000 
0.5776 
0.6431 
0.9073 
0.9828 


ll' 


1.500 


1.357 


1.483 


1.556 


1.483 


1.483 



Note: p = 0.5, True W(x) = 1 - 0.5532e _0 • 35B8x + 0.05317e" IMO1 . 

Table VI — P(Delay < x) for case M/E 2 /1 



X 


True 


A EM 


A t -., 


■Ao.1 


A* 


A*' 




0.5 

1.0 

5.0 

10.0 


0.1000 
0.1244 
0.1509 
0.3497 
0.5359 


0.1057 
0.1354 
0.1641 
0.3617 
0.5446 


0.0919 
0.1220 
0.1511 
0.3518 
0.5374 


0.1000 
0.1278 
0.1548 
0.3426 
0.5198 


0.1057 
0.1349 
0.1632 
0.3584 
0.5397 


0.1000 
0.1296 
0.1582 
0.3557 
0.5388 


w 


13.50 


13.26 


13.46 


14.33 


13.46 


13.46 



Note: p = 0.9, True W{x) = 1 - 0.9 lie- 006 ' 45 * + O.OlllOe-' 483 *. 

QUEUEING SYSTEM APPROXIMATIONS 313 



where d is the constant interarrival time. Thus, 
K(0) = B(d) 



KM) = 



J e ay dK{y) = 1 - | e ay dK(y) = 1 - K + (a) 



lu = xdK(x) = K' + (0) , 

Jo 

where a satisfies 

e ay dK(y) = K-(0) + K o (0) = I. 



This class of systems (D/G/l) is particularly important in the 
analysis of a certain class of schedules for computer systems with real 
time applications. For example, see Refs. 1 and 2 where the approxi- 
mations given here are applied to the study of these schedules. The 
service time distribution for these systems are generally discrete in 
nature. Here we look at a simple special case with 

A(t) = U(t - d) 

B{t) = p x U(t - si) + p 2 U(t - s 2 ) , 

where U(x) is the unit step function, pi + p 2 = 1 and Si<d<s 2 . From 
these, we readily obtain 

K(0)= Pl 
K-(a) = Pie a{s >- d) 
K + (a) = p 2 e a ^- d) 

Tu = K' + (0) = p 2 (s 2 - d) , 
and a.E satisfies 

K-(a) + K + (a) = e ad ( Pi e as > + p 2 e a!<2 ) = 1 . 
Thus, (7) and (9) readily yield 

Pu = Co(a E ) = l _ l ~J^_ d) (63) 

Ci(a E ) P2 (s 2 - d) 

WE,\ = = a F (8,-d) ' (b4) 

dE pi-pie aEK ' 

For simplicity, we consider the special case Si = 0, s 2 = 2d. In this 
case, the exact solution satisfies the difference equation 

Pi = piP, + i + PiPi-i, (65) 
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where P, = P(work in the system = id at an arbitrary arrival epoch), 
i.e., Pi = P(arrival delayed id). Taking d = 1, the solution to (65) is 
readily found to be 



Pi = P(Delay = i) = (1 - r)r'; r = ^. (66) 



Pi 



Hence, 





Pd 


Pi 






w 


r 


P2 




1 - r 


Pi ~P2 


and oe = ln(r 


! ). 






Now (63) and (64) yield 








C (a E ) 


I -Pi 
l-p 2 


= P2 
Pi 






d(a E ) 


P2 



(67) 
(68) 



a E Pi - p 2 



(69) 
(70) 



i.e., Ae.o gives the exact value for Pd and Ae.i gives the exacf value id. 
In particular, this implies that the approximation Ae.o is exact at grid 
points — Wa B9 (x) is exact for x = i, an integer. 

Table VII shows the resulting delay distributions for several of our 
approximations. We see here that A* seems to be the best overall 
approximation. Note that A ,i (and hence A**) result in somewhat 
poorer approximations for this case. This may indicate that A* is 
potentially more robust. 

As another comparison for the mean delays predicted by our ap- 
proximations, we have included on the table (denoted by w[K — L]) 
the results of using an approximation for the mean delay given by 

Table VII — P(Delay < /') for case D/D 2 /1 



i 


True 


Ae.O 


A E ,\ 


A* 


A ,\ 





0.1818 


0.1818 


0.0970 


0.1818 


0.1282 


1 


0.3306 


0.3306 


0.2612 


0.3178 


0.2331 


2 


0.4523 


0.4523 


0.3955 


0.4312 


0.3254 


3 


0.5519 


0.5519 


0.5054 


0.5258 


0.4066 


4 


0.6334 


0.6334 


0.5953 


0.6046 


0.4780 


5 


0.7000 


0.7000 


0.6689 


0.6704 


0.5408 


6 


0.7546 


0.7546 


0.7291 


0.7252 


0.5960 


7 


0.7992 


0.7992 


0.7784 


0.7709 


0.6446 


8 


0.8357 


0.8357 


0.8187 


0.8089 


0.6874 


9 


0.8655 


0.8655 


0.8516 


0.8407 


0.7250 


10 


0.8900 


0.8900 


0.8786 


0.8672 


0.7581 


w 


4.50 


4.08 


4.50 


4.50 


6.80 


Note: p 


= 0.9, True P( Delay - i) = 1 1 


-(tt)1 


(tt) 


, w(K - L) = 6.43. 








w. 


\ll/ 
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Kramer and Largenbach-Belz. 11 This approximation is essentially a 
heuristic extension of an approximation originally obtained by Hey- 
man 12 via diffusion techniques. The rather simple form of the approx- 
imation w(K - L) given below makes its use very appealing. 

w(K - L) = 



2(1 -p) 



where 



Y = 



2(l-p)(l-C 2 ) 2 2 

3p(C T 2 + C 2 ) 
(1 - p)(C 2 - 1) 2 



and C T , C t are, respectively, the coefficients of variation of the service 
time and the interarrival time. (Note that w(K - L) is exact for an 
M/G/l system.) 

6.4 Case:H 2 /E 3 /1 

As a last numerical example, we consider an H 2 /E 3 /l system. Spe- 
cifically, 



B(s) = 



Xi + s X2 + s 

1 



(1 + s) 3 

Again, the exact waiting time is relatively easy to find via standard 
techniques (see Ref. 10), which yield 

W(x) = 1 - Ce-"" - C 2 e- a2X - C 2 e-°* x , 
where the a, are the roots of the equation 

R(s)=A(-s)B(s) = 1, 

a denotes the complex conjugate of a, and the C,'s are the correspond- 
ing residues which are readily obtainable (see Ref. 10). 
For the approximations of the preceding sections, we find that 



P\ , P2\ ( Pi , Pi 



°« = 3 + 2 'x! + xi - x; + f 



K(0) =p t q} +Piq2 

& , x Piqi^i . P292A2 

K-(a) = -— + t— , 

Ai + a\ A2 + a.\ 
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Table VIM— P(Delay < x) for case H 2 /E 3 /1 



* 


True 


A* 


A** 


Ao,, 


A'o,, 





0.7834 


0.7897 


0.7835 


0.7835 


0.7835 


0.5 


0.8121 


0.8265 


0.8225 


0.8212 


0.8080 


1 


0.8406 


0.8569 


0.8544 


0.8523 


0.8377 


5 


0.9731 


0.9694 


0.9702 


0.9680 


0.9740 


10 


0.9977 


0.9955 


0.9959 


0.9953 


0.9978 


20 


0.9999844 


0.9999054 


0.9999223 


0.9998966 


0.9999849 


w 


0.5536 


0.5458 


0.5458 


0.5662 


0.5523 



Note: p = 0.2, True P(Delay < *) = 1.0 - 0.3260e" 04974 ' + e _12,5l [0. 10944 cos(0.3129*) 
0.08673 sin (0.3129*) ], W Aoi (x) = 1 - O^ieie -04974 * + 0.0995e _1 2i5x , w(K - L) = 0.5476. 



Table IX — P(Delay < x) for case H 2 /E 3 /1 

True A* A" A^ t A^ 

0.5867 0.5730 0.5730 0.5727 

0.6408 0.6305 0.6280 0.6104 

0.6877 0.6803 0.6752 0.6550 

0.8982 0.8996 0.8913 0.9008 

0.9749 0.9764 0.9723 0.9808 

0.9985 0.9987 0.9982 0.9993 






0.5728 


0.5 


0.6145 


1 


0.6578 


5 


0.8997 


10 


0.9806 


20 


0.9993 



40 


0.9999990 


0.99999997 


0.9999998 


0.9969756 


0.9999990 


w 


1.497 


1.475 


1.475 


1.561 


1.494 



Note: p = 0.4, True P(Delay < x) = 1.0 - 0.5176e-° 32Ml + e~' 265 T0.09038 cos(0.37! 
- 0.06644 sin (0.3793*)], W A Jx) = 1 - 0.5124e-° 3285 " + 0.0851e- ,26Sl , w(K - L) = 1. 

Table X — P(Delay < x) for case H 2 /E 3 /1 

x True A* ~A** ~A^ A ^7 



3793*) 
488. 






0.1785 


0.1906 


0.1791 


0.1791 


0.1785 


0.5 


0.2041 


0.2253 


0.2148 


0.2122 


0.2024 


1.0 


0.2325 


0.2586 


0.2489 


0.2440 


0.2314 


5.0 


0.4609 


0.4779 


0.4738 


0.4562 


0.4616 


10.0 


0.6588 


0.6632 


0.6627 


0.6398 


0.6592 


50.0 


0.9912 


0.9899 


0.9904 


0.9866 


0.9912 



100 


0.9999093 


0.9998741 


0.9998874 


0.9997827 


0.9999094 


w 


9.285 


9.230 


9.230 


9.967 


9.928 



Note: p = 0.8, True P(Delay < *) = 1.0 - 0.8517e" OO915j + e~' 32,I [0.03026 cos(0.4548*) 
- 0.02125 sin (0.4548*)], W A Jx) ;q 1 - 0.8517 - 009 ' 5 * + 0.0302e-' 32U , w(K - L) = 9.266. 

where 

Tables VIII, IX, and X compare our various approximations for this 
case with p x = 0.25, A 2 = 2A, and p = 0.2, 0.4, and 0.8. + (The approxi- 
mation denoted by A'o,] will be discussed shortly.) Here we see again 
that A* and A** yield quite good approximations over a wide range. 

Again, for comparison, we have included the approximation w(K — 



f The roots needed for the exact solution (and the approximations) were obtained by 
using a program developed by A. E. Eckberg. 
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L). Unlike the D/D2/I case, we see that w(K - L) provides a good 
approximation to the mean delay for these cases. 

VII. EXTENSIONS 

As we have seen, even with the assumption of only a single expo- 
nential for the form of W(x), we are led to a wide variety of approxi- 
mations by using the Lindley integral equation (1) to determine the 
two unknown coefficients. We consider here some extensions for the 
case where we wish to choose a more complicated form for W{x). 

7. 1 Method of moments for hyperexponentials 
In many cases, W(x) admits the expansion 

W(x) - 1 - S c,e-°' x , (71) 



1-1 



where the coefficients C, a, are not necessarily real. In such cases, an 
approximation of the form 

m 

W A (x) = 1 - 2 &e- aiX (72) 

i-l 

seems most appropriate. In general, of course, we do not know the 
form of W{x). However, if the equation 

g(-s) = A(s)B(-s) = 1 (73) 

can be solved more than one nonzero root, then it seems reasonable to 
attempt to incorporate additional roots in an approximation of the 
form (72). One can discretize (1) to obtain a set of (implicit) equations 
for the needed C; however, we show here that the method of moments 
introduced in Section III can readily be extended to obtain a set of 
linear equations for these quantities. For this purpose, it is convenient 
to consider the equation for W c {x) = 1 - W(x) corresponding to (1), 
i.e., 

L[W c ](x) = W c (x) - 1 + K(x) - J W c (x - y)dK(y) = 0. (74) 

We denote the Lth moment of (74) by n L : 

H - J x L dL[W c ](x). (75) 

Jo 



Using (74) in (75) we obtain 



= J dW c ( 

Jo 



(x) + dK(x) + W c (-y)dK(y) 



Mo 

'0 
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= j x L dW c (x) + I 
Jo Jo 



fiL = x L dW c (x) + x L dK(x) 



+ L\ x'" 1 \ W(x - y)dK(y)dx, L>0. (76) 

'0 J-m 



With the notation 

W L = | x L d\\"(x), u, - | x L dKU) 
(76) can be written 

H = wl + m+ | it"( y)dK(y) 

Hl = wi + uL 

+ L\ x L - } I W c (x-y)dK(y)dx, L > 0. (77) 



I x L dW c (x), ui =■ I 
Jo Jo 

f. 



[Note Mo = 1 - #(0) and «i = m + of the preceding sections.] 
Now using the form 

Wa(x) = S C,e°<* (78) 

/io = leads to 

2 CW - 1) = -55, (79) 

where 

K,= e°- y dK(y). 



After some algebra, ixl = for L > yields 

«, „ (UKj ^ Lu L -i ^ y L(L - 1 ) ... (L-*) _ L! 

* C, l~^" + _ ^~ + ? 7F* UL-l-k--Z 

, \ a, a, /,-i o.i a , 

= -u L . (80) 

The empty sum is taken to be zero. 

Thus, given a set of roots a„ i = 1, ■ • • , m of (73) eqs. (79) and (80) 
allow us to readily compute the desired coefficient C,-. (Note that this 
method is identical to the common techniques of "method of moments" 
frequently used for other classes of integral equations [e.g., see Ref. 
13).] We refer to the resulting approximation as 

AJ,....^. 
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Table XI — Dominant residue for case H 2 /Ei/1 



p 


True 


Cejq 


Cea 


Co.i 


C t (A ,',,,) 


0.2 

0.4 
0.8 


0.3260 
0.5176 
0.8517 


0.2103 
0.4133 
0.8094 


0.2714 
0.4841 
0.8444 


0.2164 
0.4270 
0.8209 


0.3161 
0.5124 
0.8517 



7.2 Numerical example (H 2 /E 3 /1) 

Note that if we apply the results of (71) to the H 2 /E 3 /l example of 
Section VI and include all three roots, a u a 2 , a 3 , then the resulting 
delay distribution will be exact [that is, if all roots are given, (80) is 
equivalent to the standard methods of determining the appropriate 
residues]. To see how we can use more structure to improve our 
approximations, we will use the true dominant root ai as one of our 
roots, but take a' 2 = Re(a 2 ) as another root, the resulting approxima- 
tion is shown on Tables VIII, IX, and X, where it is denoted by Ao,i. 
We see that the inclusion of this additional term results in an approx- 
imation that more closely captures the structure of the true delay 
distribution. This is perhaps best illustrated by Table XI. Here we 
have given the values for the dominant residue, i.e., the coefficient of 
the exponential with the dominant exponent, for the various approxi- 
mations. We see that of the three single exponential approximations 
A E ,\ has a C value nearest the true value. This shows why A E ,\ tends 
to have better tail behavior than the other single exponential approx- 
imation. However, we see that Ao.i resulted in an excellent approxi- 
mation to the true dominant residue, even though we did not use 
another exact root. Hence, the excellent agreement on Tables VIII to 
X for the extreme tails of the distribution. This behavior can be very 
important in studying computer systems with dedicated real-time 
applications where often criteria are specified in the 10 -5 probability 
range, i.e., the probability of delay greater than T shall be less than 
10 -5 . 

VIII. CONCLUSIONS 

The basic idea we have exploited is to choose a functional form for 
an approximation W A (x) to a true delay distribution, say, W(x), and 
use the well-known Lindley integral equation to find the undetermined 
coefficients. For the case where W A (x) is exponential— W A (x) = 1 - 
Ce~°* — we have used this technique to develop several approximations, 
some of which make use of the explicit structure of the relevant service 
and interarrival time distributions, while others require only moment 
information. 

Although not always the best choice, A* seems to provide the 
robustness that one would require of a good approximation. The 
resulting approximation for the mean delay is excellent and the result- 
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ing probability of delay quite good, although A** provides a better 
value for the latter. For predicting tails of distributions, we see that 
Ae,\ is the best of the simple exponential approximations. We have 
also seen that increasing the complexity of the forms of the waiting 
time distribution assumed in the approximation, e.g., using more than 
one exponential, can result in extremely accurate predictions of the 
tails of the delay distribution. 
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APPENDIX 

Summary of Notations and Formulas 

t = interarrival time 
A(t) = interarrival time distribution 
L.S.T. = Laplace-Stieltes Transform 
A(s) = L.S.T. of A(t) 
E = Expected value 
t = E(i) = 1/A 
t = service time 
B(t) = service time distribution 
B(s) = L.S.T. of B(t) 
t = E(t) = l/o 
P = A/a m 

U = T — t 

K(u) = distribution function of u 
K(s) = L.S.T. of K(u) = A(-s)B(s) 
u = E(u) = 1/a- 1/A 
ol — variance of u 

K_( S ) = e m dK(u) 



K + (s) = e su dK(u) 
Jo 

K(s) = K-(s) + K + (s) 
a E = positive real root of I e au dK(u) = 1 
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a, = ordered roots of characteristic equation, 
K(-a) =A(s)B(-s) = l(a, = a E ) 

u + =\ xdK(x) = K' + (0) 







K(0) = K-(0) 
an = —2u/al 
p=A(a) 
q = B(X) 
W(x) = waiting time distribution 

P D = true probability delay greater than zero 
u = true mean delay 
W c (x) = 1 - W(x) = complimentary waiting time distribution 
L[W c ](x) = Lindley integral equation (complimentary) 



y)dK(y) 



L[W c ]{x) = W c (x) - 1 + K(x) - J W c (x - 

H = \ x L dL[W c ](x) 

Jo 

Ki = K-(ai) 
Approximate waiting time distribution (single exponential): 

W A (x) = 1 - Ce~ ax 

C = approximate probability of delay 

w A = C/a = approximate mean delay . 

Approximate A E .o- 

a = a E 

n n v 1 - *<°> I ~ ^- (0) 

C ~ C ° iaE) ~l-Ua E )~l-K-(a E )- 

Approximation AeJ. 

a = a E 

Ci(a E ) u + K' + {0) 



we.x = 



a E K(0) - K-(a E ) K-(0) - K-(a E ) 
Approximation Ao,H 
a = ao.i solution of 

[(K(0) - K-(a)][l - K(0)] _ [K-(0) - £-(a)][l - £-(0)] 
a= u+[l-K-(a)] £' + (0)[l - £-<a)] 
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_ n 1-KiO) l-K-(0) 



l-K-(ao A ) 1-K-(ooa) 
Approximation A*: 

C=C* = C (a E ) (fromA £ ,o) 

_ + C* CMe) 

W m - — - = (from A E ,i) 

a OLE 

, . » Co(a E ) 

determines a * = — — a^ . 

Ci(a£) 

Approximation A**: 

C=C** = C .i (fromAo.i) 

_„ C** d(g g ) #iL 

w** = ^rr = (from^.i) 

a** as 

determines a * * = ' a^ . 
Ci(a£:) 

Approximation Ah,q' 

-2m 



a = a H = — o" 
ol 



C = C//,o = Co(a./i) = 



1 - K-(a H ) 1 - /£-(a ff ) 
Approximation Aha' 
a = a„ 

Cha CAa„) u + K' + (0) 



Wha 



a H a„ K(0) - K-(a H ) £-(0) - K-(a H ) ' 

Approximation Aht: 

a = a H 

Cht 1 

Wht = = — . 

a« oh 



Approximation Ag,o' 



a = an 



Cg.o = 



VvW 



1 + Erf( - " 



\y/2a u , 
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Approximation Aq,\\ 
a = a H 



o u 4(r) . u 



2 \o„ 

Cg.i J2tt 
w g ,\ = 



Approximation Alt: 

Cat * 1 - *<0) = 1 - ^-(0) 

U ' LT " a/.r ^(0) £-(0) " 
Approximate waiting time distribution (multiple exponentials): 

m 

W A (x) -i-2 G«r°* 

1-1 
Approximation Aq,i,. ., , m ' 

a, roots of K(-a) = 1 
al -at, i = 1, • • • , m 
d determined from 
S ^ ( UKi ^ La ^ j. V L(L-l)---(L-k) _L\\_ 

,_i \ a; a, A-i a, a. / 

where the empty sum is taken to be zero. 
Approximation A'oy. 

ai roots of K(-s) =A(s)B(-s) = 1 

a\ = a\ = a E 

a'i = Re(a 2 ) 

d determined from 

(1 - K x )Cx + (1 - K 2 )C 2 = m = 1 - K(0) 

a'i a2 
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